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space dimension. The exact reactive (fluctuating) steady state and the long-time 
behavior of the concentration of particles are explicitly computed. 
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1. Introduction 

In the last decades much effort has been devoted to the study of nonequilibrium 
systems, both theoretically (see e.g. and references therein) and experimentally 

(see e.g. 01 El HH HI an d references therein). Despite these efforts, a comprehensive 
theoretical framework is still lacking: As yet, there is no equivalent of ensemble theory 
for nonequilibrium systems. Consequently, in nonequilibrium systems one has to 
explicitly impose the dynamics in order to be able to compute observable quantities. 
In this context the stochastic dynamics is generally described by a master equation, 
and most progress in the field is made by studying paradigmatic models pQ. In this 
framework, exact solutions of simple models are scarce, but very precious, since they 
can serve as benchmarks for approximate and/or numerical schemes and shed light 
on general properties of whole classes of related models. Not surprisingly, nontrivial 
solutions are almost entirely restricted to one dimension (ID; see e.g. 0G] ) and have 
focused mainly on spatially homogeneous models. One should however mention that 
the authors of Ref. |Hj were able to obtain the exact steady-state of a one-dimensional 
coagulation-diffusion model with local particle input, while the author of Ref. [5] 
computed the non-trivial stationary concentration of particles in a one-dimensional 
system of reversible diffusion-limited coalescence in the presence of a (perfect) static 
trap. Also, recently the inhomogeneous version of the voter model, and its mapping 
to the monomer-monomer catalytic reaction have been solved in arbitrary dimensions 

nniE]. 

Actually, in many realistic situations, diffusion of particles (describing e.g. ions 
or excitons mEJE]) often takes place at interfaces of different phases (inhomogeneous 
substrates) and thus the understanding of the spatiotemporal effects of impurities, 
defects, and inhomogeneities is highly desirable. To gain some further insight 
in these issues, the analysis of simple models seems to be the most natural and 
promising approach. Motivated by these considerations, and complementary to recent 
studies on inhomogeneous reaction-diffusion systems (see e.g. [HHH]), we specifically 
investigate, through exact methods and in arbitrary dimensions, the properties of an 
inhomogeneous diffusive stochastic many-body system where particles are injected 
by a uniform source which refills empty sites of a lattice. Once on the lattice the 
particles can jump symmetrically onto their empty neighboring sites. We also consider 
that at the origin of the lattice there is trap that can remove particles from the 
system. The goal here is to carefully analyze the effect of the competition between 
the uniform source and the local trapping reaction in various dimensions and to study 
the spatiotemporal concentration profile of particles around the trap. It is therefore 
shown that the interplay of these competing (homogeneous vs local) reactions results 
in nontrivial stationary states and affects also the time behavior of the system. In 
addition, as there have been recently various experimental works d devoted to the 
study of the properties of such a model in the absence of the source, it seems that our 
study could also have some experimental implications as well. 

This work is organized as follows: In the next section we give a precise 
mathematical definition of the model under consideration, the basic equation and 
an overview of the main results. In Section 3 we set-up the formalism to compute 
the dynamical concentration of particles. In Section 4, we especially focus on the 
ID infinite chain and obtain explicitly the stationary and long-time behavior of the 
concentration. In Section 5, we focus on the two-dimensional infinite system, while 
the section 6 deals with the three-dimensional case (in the thermodynamic limit). In 
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Table 1. The three elementary processes under consideration. As explained 
in the text, the hopping and injection processes are homogeneous, whereas the 
trapping reaction is localized at site "0" . 

Reactions Processes 

A0 < — ' 0A Symmetric hopping with exclusion 

-^-> A Homogeneous source refilling empty sites 

A > Trapping of particles at site "0" 




Figure 1. Illustration of the processes occurring in the one-dimensional version 
of the model, a) Particles can jump to neighboring empty site (with rate D) and 
can be injected on empty sites by a source of intensity J. The site at the center 
of the lattice is the "trap" and the particles (indicated in black) landing on this 
site can be removed with a rate T. b) Configuration of the system after all the 
moves sketched in a) have been performed. 

Section 7, we study the global effect of the trap and the source for infinite systems 
in ID, 2D and 3D. In Appendix A, we briefly review the results already obtained 
(through other methods [1211131 H^lll5p for the diffusive model in the presence of the 
trap but in the absence of the source (and in the thermodynamic limit). In Appendix 
B, one extends the discussion to the case of periodic boundary conditions on finite 
lattices and especially focuses on the concentration of particles on a ring (ID) of finite 
size. 

2. The Model and outline of the results 

In dimensions d — 1,2 and 3, we consider a stochastic model of interacting particles on 
a hypercubic lattice of size (2L + l) d . Each site, that can at most be occupied by one 
particle, is labelled by a vector x of components — L < xi < L, i — 1, . . . , d. On this 
lattice, particles denoted A can jump (subject to the site restriction) to a neighboring 
empty site (denoted by the symbol 0) with a rate D and an uniform source refills 
the empty sites of the lattice by injecting particles with rate J. In addition to these 
homogeneous processes, we consider that there is a trap at the origin of the lattice such 
that a particle landing on this "poisoned" site can be trapped and removed from the 



Competition between a uniform source and a local trap in a stochastic system 



4 



system with a rate T. These processes are summarized in the Table 1 and illustrated 
in the Figure 1. 

The competition between the homogeneous process refilling the system and the 
local one trapping the particles is studied by investigating the properties of the 
concentration of particles. The equation of motion for such a quantity is obtained 
from the master equation and reads : 

= D A x c x (t) + J[l - c x (t)] - Tco(t)8 a , . (1) 

at 

Here denotes the discrete Laplace operator: A x c x (t) = —2dc x (t) + Ylx' Cx ' (*) > 
where the sum on right-hand side (r.h.s.) runs over the 2d nearest neighbors x' of site 
x. It should be noted that the equation (JJJ is both valid for infinite systems (L — > oo) 
and for periodic lattices of arbitrary size. In Sections 3 to 7, for the sake of technical 
simplicity, Eq. Q will be studied in the thermodynamic limit (L — > oo). The case of 
periodic boundary condition, somewhat technically more involved, is briefly addressed 
in Appendix B. 

At this point, some comments on the diffusion-like Eq. Q are in order. We 
first notice that the effect of the trap appears through the inhomogeneous term 
— Tca(t)5 Wi o which is non zero only at the origin of the lattice and is proportional to the 
concentration. In addition to the inhomogeneity, the equation is not merely diffusive 
since there is a term accounting for the injection of particles by the source. Because of 
the site restriction the system allows at most one particle at each site. Consequently 
the source term refills only empty sites and contributes with a term J[l — c x {t)\ (and 
not simply J). One can also mention that systems with site restrictions are usually 
not solvable in arbitrary dimensions: when available, most of the exact solutions 
are obtained in one-dimension (see e.g. and references therein). It is also worth 

noting that the authors of Refs. [HUH] were able to solve some spatially inhomogeneous 
versions of the coagulation-diffusion stochastic model in one dimension by using the 
so-called IPDF method (see e.g. and references therein). 

Remarkably, for the model under consideration, the equation for the 
concentration is linear (but inhomogeneous) and thus solvable. Obviously, in the 
absence of the source term, the system reduces to a purely diffusive model in the 
presence of a trap and, because the hopping process is symmetric, the volume exclusion 
does not affect the concentration of the particles [see Eq.QJ]. In this case one 
recovers the results obtained in previous studies based on non-interacting random 
walks jHJ |3 ^3 1141 ITS] , as it is briefly reviewed (for a discrete lattice) in Appendix 
A. Also, it is clear that in the absence of the trapping reaction (T = and J > 0), 
the system is eventually completely occupied by the particles and (for homogeneous 
initial conditions), the approach toward the steady state is purely exponential in all 
dimensions: c x (t) — 1 ~ e~ Jt [see Eq.Q]. In the presence of the source, the role of the 
site restriction does effectively matter. In fact, if one would relax such a constraint [and 
would have J instead of J(l — c x (t)) in Eq.JTJ], the concentration of particles would 
simply increase linearly with the time and there would be no effective competition 
with the trapping reaction. In this sense, the many-body character of the system 
under consideration is embodied in Eq.QJ by the term J(l — c x (t)): as the source 
injects particles only on empty sites, the interplay with the trap drives the system 
toward a non-trivial genuine nonequilibrium reactive (fluctuating) steady state which 
properties explicitly depend on the dimension of the system. Of course, the most 
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interesting situation arises when the homogeneous injection rate is small (but not 
zero) and is thus efficiently balanced by the trapping rate. In fact, as we show in the 
next sections, the effect of the competition between the local/homogeneous reactions 
(on infinite systems) results in a perturbation of the uniform static density profile 
Sc x (oo) = 1 — c x (oo) which reads (x = \x\) 



Sc x (oo) = < 



2K (x^/j/D) 
In (D/J) 

[Aze-*Vm/x 



; in ID 
; in 2D 
; in 3D, 



(2) 



where £ = 2 In ( ^+^l+ 4g ) < ^/J/D, K q (z) is the Bessel function of third kind 

and Ai^ are amplitudes which values are given by Eqs JEJ) and J2ZJ. It is worth 
mentioning that the ID and 3D results quoted in Eq.Q are valid for any value of the 
parameters. However, while the ID expression is the exact solution of the stationary 
difference equation the 3D result has been obtained from a suitable continuum 
reformulation of Eq.JQl. On the other hand, the above 2D expression of 5c x (oo) has 
been obtained from Eq.QJ in the limit where J is small. 

We also show that the time evolution of the concentration is affected by the 
competition between the source and the trap and depends on the dimensionality of 
system. In fact, it is found (on infinite systems) that the time-dependence is not 
purely exponential and depends on the initial state of the lattice: 

-Jt, 



-,(t) - c x (co) oc < e J (lni) 



-Jt t 



; in ID 
; in 2D 
; in 3D, 



(3) 



where the subdominant time-dependent prefactors have the exponents v\ = vz = 
3/2 and r\ = 2 for an initially fully occupied lattice and v\ = v% = 1/2, rj = 1 for other 
(homogeneous) initial conditions. All these results are derived and discussed in detail 
in the remainder of this work. 



3. General set-up 

Using the properties of the Bessel functions of first kind I n {z) = Jq ^■e zcosq cos nq, 
[TC] and extending the techniques devised and briefly presented in Ref. JO] (see also 
one derives from equation the following exact and self-consistent relationship 
obeyed by the concentration in the thermodynamic limit (L — > oo): 



c x (t) = l+5>«(0) - l)e- Jt f[[e-* Dt I x ^ yi (2Dt)] 

V 8=1 

-T f dt'c (t - t') e- n ' TT \e- 2Dt 'l Xt (2Dt')} , (4) 

where the sum on the r.h.s takes care of the initial condition and the convolution term 
accounts for the inhomogeneity. Of course, this expression is valid both in the absence 
(J = 0) and in the presence (J > 0) of the source. It is clear from this expression 
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that in the absence of the trapping reaction (i.e. if T = Q), at long time the lattice 
is fully occupied: c x (oo) = 1. Obviously, the situation is much more interesting 
and subtle when there is a competition between the trap and the source. It is also 
important to notice that the explicit determination of the density of particles still 
requires the solution of the self-consistent integral equation (@J, where the expression 
of the concentration appears both on the right (r.h.s) and left hand-side (l.h.s). This 
is achieved by using Laplace transform techniques and the integral representation of 
the modified Bessel functions. In fact, we introduce the quantity T x that is essentially 
the Laplace transform of a product of Bessel functions (multiplied by an exponential 
factor): 

2.M,J)=/ dte-W* T\{e- 2m I Xi {2Dt)} . (5) 

From the integral representation of the modified Bessel functions, one can check that 
the quantities I x {s, D, J), hereafter simply denoted T x (s), can be rewritten as the 
so-called Watson integrals (or cubic lattice Green- functions) : 

1-*"']-* {^Ys + J + 2D£ii( 1 -<x»«)' 

(6) 

where we have introduced a <i-dimensional vector q — (q\, . . . , qd). For further 
convenience we introduce the norm of the vector x, that is denoted x = \x\ = 



Si=i x i ^ 0, and the parameter u = xyj J/D. Also, as in arbitrary dimensions 
the concentration of particles depends spatially only on the distance x to the origin 
[see e.g. Eq. hereafter it is simply denoted c x {t) (instead of c x (t)). 

We now focus on the site x = 0, and specifically consider that initially the system 
is randomly occupied with an homogeneous concentration c x (0) = po of particles. 
Using the convolution theorem and the identity J2n=-oc ^n{ z ) = e z from Eq. Q 
we get the Laplace transform of the concentration of particles at the origin: 



DC 



Plugging back the result Q into the Laplace transform of Equation (0J, we obtain: 



c x {s) = dte st c x (t) = 



l+T[I (s)-I x (s)} 



1 1-po 



J 



(8) 



l + TT (s) 

Again, the expressions and © are valid for any strength J > of the source. 
Using the general properties of Laplace transforms |16| , the concentration for a system 
with initially c x (0) — po is simply related to the concentration obtained from the case 
where po = 1 and J = 0. In fact, it readily follows from (JJJ) and JHJ that 

CxWUo,.*,,,^ = CxWUo^x - (! - P^' 1 [c,(*)| J= o iPO=1 ] • (9) 

The long-time behavior of the concentration of particles is therefore obtained from the 
analysis of J7J and JSJ. This study is the scope of the the subsequent sections where 
due attention is paid to the competition between the homogeneous/local reactions and 
to the effect of the dimensionality on the spatiotemporal properties of the system. 
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4. One-dimensional concentration in the presence of the source and the 
trap 

In this section we compute explicitly the stationary concentration as well as its long- 
time behavior on an infinite chain (ID). 

In one dimension, the quantity T x (s) can be computed explicitly and one has 
[x > 0, see Ref. [TB]): 

V ; V(s + J)(s + J + 4D) \ 2^D J V ' 

It is clear from this expression that, as long as J > 0, T x (s) is a well defined 
quantity, even in the s — > limit. However, when J = 0, in the regime s — > 0, 
T x (s) diverges as s -1 / 2 . Because of this fact, we can already anticipate an exponential 
dynamics in the presence of the source and an algebraic one in absence of the latter 
(see Appendix A). In this section, we focus on the case where the source is non-zero. 

Clearly, it follows from JSJ , that the stationary concentration reads (for arbitrary 
initial conditions): 

T1 X (0) 

c - (oo) = 1 ' i + rx o( o) - 

In fact, the formula Ijlll) holds in any dimensions (when J > 0): The dimension 
d enters the expression (lllf) of c x {oo) through the value of the integrals © appearing 
in it. In one dimension, for J > 0, with (|llfl and UlUfl . we find that the stationary 
concentration approaches the value 1 exponentially with the distance to the origin: 

cjoo) = 1 ; (12) 

v ; T +V /7(7T4D)' v ; 

where one has introduced £ = 2 In ^+^+ 4£ ) -) > 0. One notices that £ < y/J/D 

and when J« 1 one has £ = ^/J/D [l - J/24D + 0{{ J/D) 2 )] . For a perfect trap, 
when T — * oo (with D and J finite), the formula l|12l) gives: c x (oo) = 1 — e~^ x and 
Co (00) = 0. These results indicate that around the origin the competition between the 
source and the trap generates a depletion zone of length £id — 2£ -1 . Interestingly, 
in the reaction-diffusion system AA «-► A it was also shown that a (perfect) trap 
also depletes the concentration profile with an effect decaying exponentially with the 
distance to the inhomogeneity [S]. Here, the situation is particularly interesting in 
the limit of a weak source, when J/T <C 1 and J/D <C 1, but J > 0. In this case, 
the competition between the localized trap and the uniform source is really effective 
and one shall distinguish two regimes [x > 0): (i) when < u = x(JfD) 1 / 2 <C 1, 
the concentration of particles increases linearly with the distance from the trap and 
with an amplitude proportional to yj : c^oo) = u; (ii) far away from the trap, when 
u is large enough, the effect of the local perturbation disappear exponentially with 
the distance to the origin: <5c x (oo) = e~ u . When J is small, the size of the depletion 
region is l\r> ~ (D/J) 1 ^ 2 . To describe further the static properties of this zone we 
can also introduce a "# s -distance" (xe 1D (oo)), which is the static distance from the 
trap to the point where the stationary concentration of A particles reaches the specific 
fraction 9 S of unity (i.e. of the stationary concentration in the bulk). Requiring that 
c x (oo) — 6 S on the l.h.s of Eq. (JT^Jl, we obtain ( v / D / 'J > 1 and y/T/J > 1): 

. , [D 2D 

x e 1D (oo) = s y — — , (13) 
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which means that xg 1D varies linearly with 8 S . 

We now compute the asymptotic approach to the stationary concentration and 
start with a system initially completely occupied (p = !)■ By noticing that 



c x {t) = 1 - f*dt'e- Jt '£- 



(t'), where C x [. . .] denotes the usual 



' T[I*(s)\ J=0 ] 
_l+T[Xo( s )| J=0 ]J 

inverse-Laplace-transform, the long-time behavior is obtained from the Laplace- 
inversion of the s — > expansion of the integrand of such an expression and one 
thus gets (Jt,—* oo): 



Co(t) - c (oo) 
and when < x < oo 

c x (t) - c x (oo) 



D 



-ji 



T Jt^Dt) 1 ' 2 ' 



-Jt-x 2 /4Dt 



jt^Dty/ 2 ' (14) 

where co(oo) and c x (oa) are given by (|12|l . Together with the results of Appendix A 
[see Eq. QA.ljl ]. we can immediately obtain the long-time behavior of the concentration 
for an initially random system, where < po < 1. In fact, according to 10 we obtain 
(0 < Pa < 1 an d t — > 00): 



Co{t) - c (00) 



2(1 ~ go) -. 
T 



(15) 



For x > and < po < 1 (with t — * 00), it is worth mentioning that there are, 
according to Eq. two contributions to the long-time behavior of c x (t) — c x (oo): 
one, the subdominant, is given by Eq. I|14fl and the other, which always dominates 
reads 5: 



c x (t) - c x (oo) 



(1 ~Po)e 



-ji 



erf 



-a; 2 /4Dt 



T 

2 



. (16) 



where erf (t) = f Q e z dz is the usual error function |lfij . Clearly, when < po < 1, 
the long-time behavior of the concentration could also have been obtained similarly 
as for Eq. Q . 

From the results l|14|) . I|15|l . (|16|l . it follows that in one-dimension the relaxation 
toward the stationary state is exponential, with a relaxation-time J -1 and with a 
subdominant power-law prefactor oc t - "^ ) which depends, as well as the amplitude, 
on the initial state of the system. In fact, v(po = 1) = | and ^(0 < po < 1) = \- In 
this sense, the long-time dynamics is said to be non-universal. Other one-dimensional 
reaction-diffusion systems displaying this kind of dynamical behavior in the presence 
of a source (but without inhomogeneities) were studied 3 . 



5. Two-dimensional concentration in the presence of the source and the 
trap 

In this section we consider the two-dimensional situation in the thermodynamic limit 
(L — * 00) and compute the concentration profile in the presence of the source. We first 

mention that in 2D and at x = 0, for J > 0, one has Io(s) = ^ ( s+ j+4g) K ( s +J+4o ) ' 

§ Of course, in the case where 1 — po X, this will require very long-time 
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where K(y) = C^ 2 de is the complete elliptic function of the first kind (see 

u y 1~ V 2 sin 2 $ 

e.g. It follows from the properties of this function that, as long as J > 0, To(s) 

is well defined for all the values of s, including s — > 0. From the previous remarks we 
infer the stationary concentration at the origin of the lattice : 

-l 



co(oo) 



r (ad 

1 + — tK 



(17) 



n{J + AD) \J + AD / 

Hereafter, we focus especially on the case where the source and the trap are in effective 
competition, i.e. when the rate of the homogeneous source is "small" compared to 
the strength of the trap and to the diffusion rate. In the limit of a weak source, i.e. 
when J/D < 1 and J/T < 1 with J > 0, we have K ([1 + J /4D]" 1 ) ~ ±ln(D/J) 
PS] and therefore c (oo) = T ^jB/j) i 1 + ° (V ln \ D / J V)}- This result shows that 
the concentration at the origin decays logarithmically with the strength of the source 
when the latter is small. We can also notice that for a perfect trap {T/D — ► oo and 
T/J —> oo), the concentration at the origin vanishes: co(oo) = 0. To compute the 
concentration of particles at sites x > 0, one has to evaluate I x (0). Following the 
same reasoning as in the previous section, we notice that, for J > 0, the stationary 
concentration of particles is given by 1)11(1 . where one has to evaluate the 2D lattice- 
Green functions at s = 0. In the physically interesting situation when the 
competition between the weak source and the trapping reaction is really effective, 
i.e. when J/D< 1 and J/T <C 1, with J > 0, the main contribution to T x (0) arises 
from the q — > expansion of the integrand of JHJ), and one thus obtains (for x large 
enough) : 

Ix{0) -> L L W J + Dq 2 = 2wD \ X ^) ' (18) 

where Ko(z) is the Bessel function of third kind ^B] which behavior for small and 

large arguments is respectively K (z) ^— ^> ln(l/z) and K (z) 'ii^ 1 (jij2z) x l' 2 e~ z . 

Using these properties in the expression (fTT]) together with the fact that X o (0) J ^J 1 
Iu(D/J)/AttD, one obtains the two-dimensional stationary concentration in presence 
of a weak source, which reads (for x > large enough): 



2K< 



(x^JjD 



c * (oo) = 1 — inWT) ' (19) 

We notice that this expression is independent of the strength T of the trap. Again, 
to discuss the perturbation due to the trap one has to distinguish two regimes 
(x > 0): (i) when < u = x(JfD) 1 / 2 <C 1, the concentration of particles increases 
logarithmically with the distance to the trap, with an amplitude proportional to 
1/lnJ^ 1 : c x (oo) ~ \nx 2 / ln (D j J); (ii) far away from the trap, when u is large 
enough, the effect of the local perturbation vanishes exponentially with the distance 
to the origin: Jc^oo) ~ ^Jir/2u e~ u j ln (D /J) . These results indicate that around 
the origin the competition between the (weak) source and the trap generates again 
a depletion zone (a circle) of radius I2D ~ (D/J) 1 / 2 . However the properties of the 
concentration profile are significantly different from those of the one-dimensional case. 
In fact, because of the subdominant y // ir/2u prefactor, outside of the depletion region, 
when x ^> £ 2 d, the effect of the perturbation decays faster with the distance to the trap 
than in ID (where it decays merely as e~ u ). The properties of the depletion region can 
be investigated further by computing the previously introduced static "^-distance" , 
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xg 2D (oo), which is the distance from the origin where the stationary concentration is 
the fraction 8 S , from i|19|) we obtain (J <C D and J ^T): 

This result shows that in 2D the xg 2D (oo) distance is a nonlinear function of Sl in 
contrast to what happens in ID Ijl3(l and it is smaller than in the one-dimensional 
case: xg 1D (oo) > xg 2D (oo). From (|12fl and l|19|) . we can also compare the effect of the 
trap at the edge of the one and two-dimensional depletion zones in the limit of a weak 
source (i.e. for x ~ {D/J) 1 / 2 > 1) by computing gj-4g$J^ ~ " 1/2 In J -1 > 1. This 
means that the difference of the effect of the perturbation due to the trap between the 
one and two-dimensional cases increases logarithmically with J . 

To compute the long-time behavior of the concentration within the depletion zone, 
we first consider the case where initially po = 1. In this situation, from Eq. (JSJ, we 

S — <§C 1 2 

obtain c x (s) ' — > — - ln ^+j)/d • After Laplace- inversion, we obtain: 

e -Jt 

c x (t) - c^oo) ~ c (t) - c (oo) ~ Jt 2 . (21) 

For homogeneous, yet random, initial conditions with < po < 1, the main 
contribution to the long-time approach toward the stationarity, within the depletion 
zone, is obtained from (|A.4|) . ljA.51) according ® : 



AttD e~ Jt 

eo(t) - c (oo) ~ -(1 - po) t (22) 
and for x ^ large enough, we have (0 < po < 1): 

g— Jt ^ rj,2 

c x (t) - c x (oo) ~ -(1 - po) ^ . (23) 

The long-time behaviors l|21|) . (1221) and (12;it show that the approach toward the 
stationary concentration is again exponential, with a relaxation time J -1 ^> 1. Since 
the dynamical prefactors and the amplitude depend on p , the long-time dynamics is 
sensitive to the initial state of the system. In contrast to the one-dimensional case 
these subdominant prefactors are not power-laws but logarithmic terms, which means 
that in 2D the approach toward the steady-state is slower than in the one-dimensional 
case [see lfT5)l and lfTH|l ]. 



6. Three-dimensional concentration in the presence of the source and the 
trap 

In the thermodynamic limit (L — > oo), we now consider the three-dimensional 
situation. In this case the integrals T x (s) [see Eq.©] are always well defined, both 
when s — > and/or J — > (i.e. limj^o^x(O) < oo). As a consequence, for all 
the values of J > the concentration profile evolves toward a reactive stationary 
state (with an infinite number of particles). This is related to the fact that in 
dimensions d > 3 random-walks are transient ^ I17j and therefore particles have 
a finite probability to be never trapped. 

It turns out that the most fruitful approach to deal with the 3D case is to adopt 
a continuum formulation suitably combined with some exact discrete results. In 
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the continuum limit, one can indeed take advantage of the fact that the continuum 
differential equation 

(A - k 2 ) N(x) = F(x), (24) 

where A is the continuum three-dimensional Laplacian and F(x) a known function, 
can be solved by the usual Green's function techniques : one has (see, e.g., [181 119| 1 

N(x) — — J d 3 x' 6 \ x _ x ,\ F{ x ')- As explained above, in the presence of a source 
J > 0, the stationary concentration of particles is given by Eq. (|llf> and an explicit 
expression of the concentration (on the discrete lattice) requires the computation of 
the 3D integrals T X (Q), with J > and x > 0. This is a hard task, and recently Delves 
and Joyce [201 have been able to calculate various properties of the simple cubic lattice 
Green functions ©, at x = and for arbitrary values of J. In particular, it follows 
from Eq.(6.12) of Reference [20]: 



_ 2(J + 6D)-y/J(J + 12D) 
0{> J 2 + 4SD 2 + Y1.JD 



1 3 V 2 
8'8 5l;r 



(25) 



where z = [j^r) , v = ^~5z-(^-z)VT=^] 2 and 2Fl ^ b . c - d ) denotes Gauss 

hypergeometric function |T§|. It follows from J7J), and (|25|l that the 3D stationary 
concentration of particles at the origin reads (for < po < 1): 

co(oo) - rpkwr (26) 

Again, one can check that for a perfect trap (T — * oo, with D and J finite), co(oo) = 0. 
We now turn to the study of the 3D stationary concentration for x > and J > 0. 
To do this, we adopt the abovementioned reformulation which has already proven to 
be fruitful in some cases (see ^3^] an d also Appendix A). One therefore considers 
N(x) = c(x,oo) — 1, where c(x, oo) is the solution of the stationary equations 
in the continuum limit. In so doing, one sees that N(x) obeys the equation l|24|) . 
with n = \J J/D and F(x) = j> i+tt (o) 3 3 (x), where the three-dimensional Dirac 
delta function 8 3 (x) has been introduced. From the solution of Eq. l|24|l . one infers the 
following stationary concentration of particles (for x > 0) : 



; (oo) — > c(a;,oo) = 1 



47rD[l+rX o (0)] x 

(27) 

Eq. I|27() . together with l|25|l and l|26|l . gives an explicit expression of the 3D 
concentration of particles. In fact, the result l|27|l is the exact continuum 3D stationary 
concentration of particles, valid for any values of J > and is in excellent agreement 
with the stationary discrete solution of 

The result l|27|l shows that in 3D the effect of the trap is much "softer" than 
in low dimensions: as the spatial perturbation of the static concentration due to the 
trap is oc e~ x ^J J I D /x, there is only a local deviation of the concentration from the 
value one rather than a real depletion, contrary to what happens in low dimensions. 
In fact, the properties of the 3D concentration profile significantly differ from those 
obtained in ID and 2D: (i) for a weak source, i.e. J/D <C 1 and J/T <C 1, the 
concentration approaches the value 1 as Sc x (oo) oc x , (ii) while, when u 3> 1 the 
effect of the perturbation is not merely exponential with the distance to the trap: 
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dc x (oo) oc e u /x. To compare the effect of the trap in 3D with the one and two- 
dimensional cases in the limit of a weak source, using (|12fl . I|19|) an d l|27|) . one can 
compute the ratio r 13 = ^ x »i,j/D<^i,j/T<^i^0p^ between Sc x (oo) in ID and 

3D and r23 = linii>i,j/fl<i,j/T«i ^"(oojjap between the two and three-dimensional 
cases, and one finds: ri3 ~ oc x 3> 1. These results confirm that the effect of the 
trap on the stationary concentration profile is much more important in low dimensions, 
where the random walks are recurrent, than in 3D where they are transient. 

As previously, to determine rigorously the approach toward the stationary 
concentration (|27|l one should analyze the s — > behavior of iJBJl, which is difficult. 
However, it follows from Eq.QJ that for an initially completely occupied lattice 
(po = 1) we have: c x (t) — c x (oo) ~ e~ J *t~ 3 / 2 . For < po < 1, it still follows 
from © and from ©, that c x (t) — c x (oo) ~ e~ Jt t~ 1 / 2 . We therefore concludes that 
also in 3D the perturbed long-time dynamics is sensitive to the initial state of the 
system. 



7. Global effect of the trap and the source 

In the previous sections the effects of trapping and injection reactions have been 
studied via the computation of concentration at each site of the lattice. Another 
question which is also relevant is related to the global effect of the source and the trap. 
As in the absence of the trap the system is fully occupied at t — > oo, it is interesting 
to compute the total number B(oo) of particles having been removed from the lattice 
by the trap and not having been replaced by the source. To address this issue, and 
thus to gain some insight on the global effect of the competition between the trapping 
and injection reactions, one computes: B(t) = X^=-oo ■ ■ • X^=-oo(l — c x(t))i which 
provides the number of particles that have been removed from the system (without 
having been replaced) up to a given time t. It follows from Eq. (J2J: 



B(t) = Te 



-Jt 



dr e' /T co(r). 



(28) 



Without help of any specific knowledge of co(t), as one always has < c x (t) < 1, 
we readily obtain the following bounds: < B(t) < -y(l — e~' 7 *), and thus 
< B(oo) < %. It therefore follows from 128|) . that in the long-time regime: 

(29) 



T 



B{t) = -co(oo)[l 



-jt-i 



where co(oo) has been computed in the previous sections. This shows that B(t) reaches 
the stationary value B(oo) = yCo(oo) exponentially fast, with a relaxation-time J -1 . 
Let us now focus on the study of B(oo). With (T2J|, (T7Jl, and we obtain: 



B(oo) 



T 



VI T+ ^J(J + 4D) 



T 
1 



T 



K 



n(J + 2D) \J + 4D 



T 



AD 



in ID 



2D (30) 



in 3D. 



I J[l + T2o(0)] 

It is therefore clear from (|30|l that the competition between the source and the trap 
also depends globally, in a non-trivial manner, on the spatial dimension of the system. 
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In particular, in the limit of a weak source (J/D <C 1, J/T <C 1), B(oo) diverges as 
J -1 / 2 in ID, as (Jin J) -1 in 2D and as J -1 in 3D. As a consequence, we can compute 
the way in which the intensity of the source should vanish to allow the trap to remove 
a number of particles comparable to the size of the system (L , L — > oo): it turns out 
that in ID one should have J/D ~ L~ 2 , in 2D J/D ~ (L 2 lnL)" 1 and J/T ~ L~ 3 . 

8. Discussion and Conclusion 

In this work, one has studied, in dimensions d = 1,2 and 3, mainly in the 
thermodynamic limit, the kinetics of a simple stochastic many-body system where 
particles, subject to volume exclusion, diffuse and interact with a localized trap in 
the presence of an external uniform source refilling the empty sites of the lattice. It 
has been shown that the competition between the homogeneous and local reactions 
always gives rise to non-trivial reactive (fluctuating) and genuine nonequilibrium 
steady states. 

In fact, in ID, an explicit expression for the stationary concentration profile of 
particles on an infinite chain has been obtained and it has been shown that the density 
of particles deviates from the value one exponentially with the distance to the trap 
(as oc e~~^ x ), while in 2D, in the thermodynamic limit, the concentration has been 
computed in the limit of a weak source (compared to the other processes), where the 
interplay between the source and the trap is particularly effective. Interestingly, in ID 
the steady-state displays the same qualitative behavior as the one obtained in the one- 
dimensional AA <-> A reaction-diffusion system in the presence of a perfect trap HJ. 
From these results it was inferred that in low dimensions the competition between the 
trap and the source generates a depletion zone in the concentration profile. Both in ID 
and 2D these regions are isotropic around the trap and their "radii" are proportional 
to the inverse of the square root of the injection rate (when the latter is small). 
We have also quantitatively analyzed the differences in these static profiles: while in 
ID the concentration within the depletion zone increases linearly, in 2D the spatial 
perturbation is less important and indeed depends logarithmically on the distance to 
the trap. In three dimensions, in the thermodynamic limit, it has been shown, using 
a suitable continuum limit reformulation, that the effect of the perturbation due to 

the trap results only in a deviation from the value one by a term cx e~ x ^ J ^ D /x. In 
this case there is no depletion in the concentration profile and the trapping reaction is 
much less effective than in low dimensions: it just generates a short-range deviation 
from the value one. To explain these different behaviors, one has to keep in mind 
the fact that random walks are recurrent in one and two dimensions, while they are 
transient in three dimensions (and above). Hence, in low dimensions, the particles 
on the lattice are doomed to visit the origin and thus might be trapped. This is no 
longer the case in 3D, where the particles have a finite probability to never visit the 
trapping site. The long-time approach of the concentration toward the reacting steady 
state has also been computed and it was shown that the inhomogeneous dynamics is 
non universal and no longer purely exponential: there are nontrivial subdominant 
prefactors which depend on the initial state of the system. The latter are logarithmic 
functions of the time in 2D and power-laws in ID and 3D. 

In summary, the results of this work show that already a single inhomogeneity 
may deeply, and non-trivially, affect in all dimensions the spatiotemporal properties of 
a diffusive stochastic many-body system characterized by the competition between an 
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uniform source and a local trap. As the ID version of the system under consideration 
displays the same kind of behavior as another reaction-diffusion model perturbed by 
a trap, it might be that the results obtained for this simple system could also apply 
more generally to a class of systems where uniform and local reactions compete. In 
particular, it would be interesting to determine if models like the reaction-diffusion 
process AA «-> A in the presence of a trap would behave in the same manner as the 
model studied here also in two and three dimensions (and not only in ID). 
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Appendix A. Diffusive stochastic system in the presence of a trap and in 
the absence of the source 

In the absence of the source, as the hopping process is symmetric and as there is 
no difference in between exchanging two particles or leaving them in place, the site 
restriction plays no role on the physical properties of the concentration of particles. 
Therefore, when J = 0, the concentration of particles in the system can be studied 
via a model of nonintcracting particles (simple random walkers) obeying a diffusive 
equation supplemented by a radiative term |18j. Such a study has been carried out 
in the references [121 1131 ITU I15| . Hereafter, for the sake of completeness one briefly 
outlines how these results can be derived on a discrete lattice (in the thermodynamic 
limit) by using the techniques described in this work. In addition, as original results, 
the time necessary for trapping a macroscopic number of particles (i.e comparable to 
the system size) is also computed. 



Appendix A.l. The one- dimensional case 

Here, in the limit L — > oo, we consider the one-dimensional situation without the 
source (J = 0) and in this case the system reduces to a so-called "symmetric exclusion 
process" (see 0]E)) perturbed by a trap. Because of the sole presence of the trap, and 
as particles perform basically random walks, the concentration is clearly doomed to 
vanish at long-time (t — > oo). Here, one outlines the steps allowing to recover, starting 
from a discrete formulation and using other methods, the results already reported in 
Refs. ^JEHEIE]- To address this point one analyzes Eq. {jJJ m the long-time 
regime, where s — > 0. It follows from (|10|l that Iq(s) ^> g L , and, for s — > 

with X\fs kept fixed, we have I x {s) ^5 e . Inserting these expressions into the 

general formula (00) one obtains cn(s) ^> — 5= and c x (s) s ^ B° I \ — £ — 

b W W Ty/7/D w s V i+Hf- 

Therefore, after Laplace-inversion |16j . one gets the following long-time behavior of 
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the concentration (Dt — > oo and Tt — > oo) : Co(t) = \ and more generally 



c x (t) = Po erf 



2VDtJ 



) + p e&(-+^/2) erfc 




PO 



erf 



2VDt 



-x 2 /4Dt 



T 

2 



(A.l) 



where erfc(i) = 1 — erf(f) is the complementary error function [T2|. It is worth to 
notice that this result coincides with Eq. (4) of Reference [S] (see also [121 1131 IT3| L 
The expression (|A.1|) is particularly appealing in the scaling regime where x — * oo 
and t — > oo, with cc 2 /£ kept finite. In fact it is readily checked that for x/VDt — > 0, 
one has c x (i) — > 0, whereas for large time and at sites very distant from the trap, 
i.e. x 3> v 7 -Dt 3> 1, Cj;(t) — > po- This means that there is an isotropic depletion zone 
around the trap. The symmetric fronts of this region propagate as i 1 / 2 . Inside the 
depletion area, the concentration vanishes algebraically with the time and from (|A.1|I 
one gets : 



c x {i) 



(A.2) 



sJ-KDt 

This result shows that within the depletion zone the concentration of particles vanishes 
as i -1 / 2 . Let us also recall that a quantity which is theoretically relevant |12II14| and, 
more importantly, which has been measured experimentally [H] is the so-called "0- 
distance" (xg). This is the distance from the trap to the point where the concentration 
of A particles reaches the specific fraction 9 (0 < 9 < 1) of its value in the bulk 
Requiring that c x (t) = 9 po on the l.h.s of Eq. (|A.2|I . one obtains: 

, 2D 

x g {t) = 9V7Dt - — . (A.3) 

In recent experiments (and numerical simulations), the depletion zone was indeed 
found to grow as the square root of the time (>; in agreement with the theoretical 
predictions (|A.1|I . I|A.2|> and (|A.3|) [Hj^]. From Eq.(0}, we can also notice that the 
total number of particles removed from the system by the trap grows as t 1 ^ 2 in ID, 
which means that the time T necessary to eliminate a number of particles comparable 
to the system size is T ~ L 2 , with L — > oo. 



Appendix A.2. The two-dimensional case 

We now turn to the two-dimensional situation in the absence of the source, i.e. J = 
and in the limit where L — > oo. As there is no external injection of particles and 
random walks are known to be still recurrent in 2D the final concentration of particles 
vanishes, like in ID. As in the one-dimensional case, here the issue is mainly technical 
as it aims at showing how to reproduce from a discrete formulation the spatiotemporal 
structure of the concentration profile, already obtained in Refs |12M13llT4lll5| . 

Using the fact that the complete elliptic functions and T x (s) diverge 
logarithmically as K ^ 1+s y 4 j j ^ -Ins; X (s) ^5 In (s/D), 16 , and 

(x » 1) l x {s) ^5 ^K (x^/I/D), we obtain: cq(s) and 
c x (s) — po - hWs/m • Therefore, after Laplace-inverting these expressions, one 
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gets the following long-time behavior of the concentration (Dt — * oo): 

and, Dt,Tt > x 2 > 1, 

/ \ In a; 2 . 

^ = «h^' (A ' 5) 
These results show that in two dimensions, and in the absence of the source, the 
concentration at sites x <C i 1 / 2 still vanishes at long-time. More precisely, there is 
an isotropic depletion zone of time-dependent radius £(t) ~ t 1 / 2 around the origin. 
Within this region, the concentration of particles decays logarithmically with the time 
[oc l/ln(t)], i.e. much slower than in ID [see (|A.2(1 and (| A. 5|> ] . As a consequence, at 
variance with the ID situation, because of the logarithmic long-time dependence of 
the concentration, the two-dimensional dynamics does not exhibit scaling. This point 
is illustrated further by the computation of the previously introduced "^-distance" , 
xg(t), which has also very recently been measured experimentally jZj. In fact, plugging 
c(t) = 6 p on the l.h.s of (|A.5jl . one obtains (Dt > 1, Tt > 1): 

x e (t)~(Di) a \ a 2 = 6/2. (A.6) 

This result, which is in agreement with the theoretical results obtained for a continuum 
model 1 15) and with recent experimental and numerical results [Jj, shows that the 
^-distance at which the concentration of particle is 9 po evolves as a power-law 
characterized by the exponent ct2 — 9/2. This result should be contrasted with the 
result l|A.3|) . where the exponent is simply 1/2, independently of the value of 9. Also, 
it follows from Q that the total number of particles removed from the system by the 
trap grows as tj In t in 2D, which implies that the time T necessary to trap and remove 
a number of particles proportional to system size is T ~ L 2 lnL, with L — > oo. 

Appendix A. 3. The three-dimensional case 

The three-dimensional case, which has been less studied ^3], displays results which 
arc qualitatively different from the low-dimensional cases already analyzed. Here we 
consider the thermodynamic limit L — ► oo. 

As discussed above, in three dimensions, even in the absence of the source, i.e. 
when J = 0, the quantity I x (s = 0) is well defined. It follows from Eqs. JHJ and © 
that the stationary concentration of particles in the absence of source reads: 

TJ x (0) 



Cx(oo) = po 



1 



(A.7) 



1+T2o(0)_ 

To make this expression explicit, we take advantage of the recent results obtained 
by Glasser and Boersma [22J , who have been able to give the exact expression of 
T-x (0). Following their work, we introduce the quantity go = "g|^r r 2 (^)r 2 (|i) = 
0.505462..., where T(z) is the usual Euler's Gamma function [TJ]. According to 
Reference |22) . we also introduce a triplet (u x , v x , w x ) of rational numbers, depending 
on x, and whose value is given [for x = (1,0,0) to x = (5,5,5)] in the Table 2 of 
Reference [221 ■ With these notations, and with help of the results obtained in |22| . 
the expression l|A.7j) explicitly reads: 



.(oo) = p a 



T (u x g 2 + w x g + |f ) 



5 o (2D + T go) 



(A.. 
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We remark that for a perfect trap, i.e. when T — > oo, 1A.8|) reduces to an expression 



that is independent of D: c x (oo) = p Q 



1 ^ 



It is useful for the 



So 

sequel to have explicitly the expression of the concentration at the origin, where 
(u ,v ,w ) = (1,0,0) E2. With ifOJ) . we obtain: 

(A.9) 

Clearly, for a perfect trap (T — > oo and Z? finite), one has co(oo) = 0. As only a 
finite number of triplets (u x , v x , w x ) are given in [22] . and as their expression and 
computation are non-trivial, to gain some insight on the stationary concentration in 
a simple manner, it is fruitful to adopt a continuum reformulation, as already devised 
in References ^3 ^J. One therefore considers the quantity N(x) = c(x;oo) — p a , 
where c(x; oo) is the stationary concentration in the continuum limit. According to 
the continuum limit of the stationary version of Eq Q), N(x) obeys the differential 
equation i|24|) . with n = and F(x) = ^co(oo)S 3 (x). In this case, the equation 
<|24[1 for N(x) is reduced to the problem of determining the "electrostatic potential" 
due to the "charge" Tcq(oo)/D located at the origin. The "charge" at the origin 
is therefore got from co(oo) (fA.Qfl and one obtains the following expression of the 
stationary concentration in the continuum limit (see [101 111] ): 

r r 



c{x, oo) = po 



1 - 



{x > 0). (A.10) 



2tt(2£> + Tg ) 

It has been checked that this continuum expression is in excellent agreement with the 
exact discrete result (|A.8|l and even "coincides" with the latter for x 3> 1. It follows 
from this discussion that in 3D, for x > and J = 0, the stationary concentration 
is a smooth and isotropic function that evolves (for x > 0) as the inverse of the 
distance to the origin [the value of the concentration at the origin is given by (|A.9[l ] : 
c x (oo) c(a;,oo) = po[l — A(D,T)/x], where the amplitude A(D,T) = 2 -n{2D+Tg ) 
depends explicitly on the strength of the trap. For a "perfect" trap, we have the same 
behavior but with an amplitude that reads A(D,T — > oo) = l/(27rgo)- At sites very 
far away from the trap, i.e x — > oo, the stationary concentration is the initial one: 
Ca;(oo) w c(cc,oo) x ^2? p _ Also, from Eq. one can show that the total number of 
particles removed from the system by the trap grows oc t in 3D. 



Appendix B. Concentration in the presence of the source and the trap on 
periodic finite size lattices 

So far, this work has focused on infinite systems. However, it is well known that 
boundary conditions may play a crucial role in nonequilibrium statistical systems. For 
instance, even for simple one-dimensional models such as the asymmetric exclusion 
process (ASEP) [see e.g. and references therein] open boundaries induce phase 
transitions. 

As the influence of boundary conditions is worth being dealt with care, in this 
appendix one briefly reformulates the problem under consideration for the case of 
periodic boundary conditions on finite size lattices. 

To emphasize the periodicity of the system, one slightly changes the notations 
and consider that the N = 2L + 1 sites are labelled by a vector x of components 
< Xi < N— 1. Again, one considers that a trap of strength T is located at the origin 
of the lattice (site 0). 
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As it was mentioned, the equation of motion for the concentration of particles 
on periodic lattices is still given by Eq.QJ. In the absence of the trap (i.e. T = 
0), the solution of Eq.Q can be obtained following Glauber 123 and thus reads: 

c x (t)\ T =o = 1 + E y (c y (0) - lK Jt n!i [e- 2Z?t E^ = _oo^-y I +n^(2^)] J where 
the summation over ni — — oo, . . . , +oo, i = 1, . . . , d accounts for the AT— periodicity 
in all the directions i. Following the same steps as in Section 3, one can therefore derive 
the following integral relationship obeyed by the concentration on periodic lattices of 
size N: 



Cx (t) = l + ^(c y (0)-l)e- Jt n 

y i=l 
r-t d 

- T / dt'c {t-t')e- Jt 'Y\ 

J ° i=l 



-2Dt 



^ Ixi-yt+mN^Dt) 
m——oo 



-2Dt 



nv — — oo 



+riiN 



(2Dt') 



(B.l) 



Again, using the convolution theorem and the properties of Laplace transform, one 
obtains the expression for the Laplace transform of the concentration of particles and 
then proceeds as in Sections 4, 5 and 6 to obtain its explicit spatiotemporal expression 
in ID, 2D and 3D, respectively. 

Appendix B.l. Concentration of particles in the presence of the source and the trap 
on a finite ring 

For the sake of concreteness and simplicity, here one illustrates how to compute 
the concentration of particles on a ring of (finite) size N = 2L + 1, with an initial 
homogeneous concentration c x (0) = po- 

Proceeding as in Section 3, it follows directly from Ea. ljB.ip that the Laplace 
transform of the concentration of particles reads (0 < x < N — 1) 

1 +TJ2n=-oo[ I nN{s) -X x+nN (s)] 



C x (s) = 



l + TY,n=-ooZnN{s) 



1 1-po 

s S + J 



As a consequence, the stationary concentration is given by (0 < x < N — 1) 
To make use of Ea. (|10|l . one writes 

oo oo 

lx+nN{s) — 2_j [Zx+nN(s) + 1(n+l)N-x(s)] ■ 



(B.2) 



(B.3) 



(B.4) 



71 — — - OO 



Using llB~4|l and Eq.JTDJl, recalling that £ = 2 In 

= -5(N-x) 



2V5 



one obtains 



Eoo 
n— - 



qq -Lx-\-nN 



(0) 



J+4D) \ 1-e 



Plugging this result into Eq. i|B.3 



one obtains the stationary concentration of particles on a periodic lattice of finite size 



c x (oo) = 1 



r ( e -e* + e -«( w — )) 



T+ V /J(J + 4D)+ T - ^JJ{J + AD) 



(B.5) 
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Plugging x — and x — N into Eq. i|B.5|l one obtains the same expression , as 
required by the periodicity (the sites and N = 2L + 1 now coincide). Clearly, we see 
that for a finite ring, there are differences with respect to the infinite chain, discussed 
in Section 4. For instance, when the strength of the source is weak (i.e. J/T < 1 
and J/D <C 1) and < x(JfD) 1 / 2 <C 1, the stationary concentration now varies as 
c x (oo) = x £ tanh (£vV/2), i.e. the slope of the stationary concentration profile in the 
vicinity of the trap increases with the size ./V of the lattice. In the thermodynamic 
limit, where N — > oo, it is easy to check that the expression l|B.3(l simplifies and, as 
expected, one recovers the expression (I12[) obtained for an infinite chain. 

To determine the long-time behavior of the concentration on a finite ring, one 
first makes use of Ecis. ljB.4D and l |10| l. and thus finds (for J = 0): 



x+nN 



001 



7=0 



i 



, 2nJv(s)|j=0 



1 



OO 

— =y 

fe=0 



("1)* 



-(x+kN)y/s/D , e -(N{k+l}-x)y/7/l 



(B, 



T 



Following the same reasoning as in Section 4, it follows from Eq. 1|B.2[1 [see also Eq.©] 
that for an initially completely occupied system (/?o = 1) the concentration of particles 



reads c x (t) = 1 — J dt'e Jt ' C 



ijv(s)| .1=0 



(t'). Using this relationship 



i+r^]~ = _ oo i„ 1 v(s)|j=o 
together with Ea. (|B.6|l . one finds the long-time behavior (0 < x < oo and t — > oo) 



Cx{t) ~ C x (o6) 



2JWnDt 



k=0 



x (x + kN)e~ ( - x+kN '> 2 / iDt + {N(k + l)-x} e W fc+1 >- 



2 /iDt 



(B.7) 



Following the same steps as in Section 4 and with help of l|B.6jl . one can also obtain 
the long-time behavior (t — > oo) of the concentration when < po < 1. According to 
Eq. ljB.2|) that there are two contributions to c x (t) — c x (oo): the subdominant one is 
given by the r.h.s. of Ea. (|B.7|l . while the dominant contribution reads 

-(l-po)e- 7 * 



c x (t) - c x (oo 

(-1 

(l-po)e- Jt ° 



k=0 



erfc 



fkN 



V 2^[Dt 



erfc 



exp 



^ N(l -t -A-) .r 



2v / -Dt 



4L>t 



exp 



(iV[l + fc] - a;) 2 



4L>i 



(B 



As for an infinite chain, it follows from Eos. (jB.7|) and IjB.Sjl that also on a (finite) 
periodic lattice the concentration approaches its stationary value exponentially fast 
with a subdominant power-law prefactor which depends on the initial state of the 
system: c x (t) — c x (oo) ~ e~ Jt t~ 3 / 2 starting from an intially completely occupied 
lattice, while c x (t) — c x (oo) ~ e~ Jt t'~ 1 ^ 2 when the initial concentration is < po < 1. 
As a difference with the infinite chain (Section 4), one notices that for a finite ring 
there are an infinite number of subdominant terms in Eqs (jB.7|l and (|B.8fl . In fact, 
one can also verify that in the thermodynamic limit (N — > oo), the expressions (|B.7fl 
and (|B.8|) simplify considerably since only the term with k = contributes. In that 
case, as expected, one recovers the results (|14[1 and <|16|) obtained for an infinite chain. 
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Clearly, one could also consider the effect of periodic boundary conditions in 2D 
(the lattice would be a torus) and 3D. However, as already illustrated by the one- 
dimensional situation, the mathematical treatment becomes somewhat involved and 
the final results deviate from those obtained in Sections 4-7 only for finite size systems. 
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